Perturbation-based moment equation approach for flow in heterogeneous porous media: applicability range and analysis of high-order terms
نویسندگان
چکیده
We present detailed comparisons between high-resolution Monte Carlo simulation (MCS) and low-order numerical solutions of stochastic moment equations (SMEs) for the first and second statistical moments of pressure. The objective is to quantify the difference between the predictions obtained from MCS and SME. Natural formations with high permeability variability and large spatial correlation scales are of special interest for underground resources (e.g. oil and water). Consequently, we focus on such formations. We investigated fields with variance of log-permeability, r2Y , from 0.1 to 3.0 and correlation scales (normalized by domain length) of 0.05 to 0.5. In order to avoid issues related to statistical convergence and resolution level, we used 9000 highly resolved realizations of permeability for MCS. We derive exact discrete forms of the statistical moment equations. Formulations based on equations written explicitly in terms of permeability (K-based) and log-transformed permeability (Y -based) are considered. The discrete forms are applicable to systems of arbitrary variance and correlation scales. However, equations governing a particular statistical moment depend on higher moments. Thus, while the moment equations are exact, they are not closed. In particular, the discrete form of the second moment of pressure includes two triplet terms that involve log-permeability (or permeability) and pressure. We combined MCS computations with full discrete SME equations to quantify the importance of the various terms that make up the moment equations. We show that second-moment solutions obtained using a loworder Y -based SME formulation are significantly better than those from K-based formulations, especially when r2Y > 1. As a result, Y -based formulations are preferred. The two triplet terms are complex functions of the variance level and correlation length. The importance (contribution) of these triplet terms increases dramatically as r2Y increases above one. We also show that one of the triplet terms is much more important than the other. When comparing K-based MCS with Y -based SME, model differences must be taken into consideration. These differences (model errors) are due to embedded assumptions and differences in implementing the discrete forms of the equations. 2003 Elsevier Science B.V. All rights reserved.
منابع مشابه
A Comparative Study on Uncertainty Quantification for Flow in Randomly Heterogeneous Media Using Monte Carlo Simulations and Conventional and KL-Based Moment-Equation Approaches
Geological formations are ubiquitously heterogeneous, and the equations that govern flow and transport in such formations can be treated as stochastic partial differential equations. The Monte Carlo method is a straightforward approach for simulating flow in heterogeneous porous media; an alternative based on the moment-equation approach has been developed in the last two decades to reduce the ...
متن کاملStochastic analysis of transient saturated flow through heterogeneous fractured porous media: A double-permeability approach
In this paper we use a double-porosity/double-permeability approach to study transient, saturated flow in heterogeneous, fractured porous media. The spatial variabilities in both the fracture and matrix continua motivate us to treat some of the fracture and matrix medium properties as stochastic processes. Hence flow through heterogeneous, fractured porous media is also amenable to stochastic a...
متن کاملAn efficient, high-order perturbation approach for flow in random porous media via Karhunen–Lo eve and polynomial expansions
In this study, we attempt to obtain higher-order solutions of the means and (co)variances of hydraulic head for saturated flow in randomly heterogeneous porous media on the basis of the combination of Karhunen–Lo eve decomposition, polynomial expansion, and perturbation methods. We first decompose the log hydraulic conductivity Y 1⁄4 lnKs as an infinite series on the basis of a set of orthogona...
متن کاملA “v2-f Based” Macroscopic K-Ε Model for Turbulent Flow through Porous Media
In this paper a new macroscopic k-ε model is developed and validated for turbulent flow through porous media for a wide range of porosities. The morphology of porous media is simulated by a periodic array of square cylinders. In the first step, calculations based on microscopic v2 − f model are conducted using a Galerkin/Least-Squares finite element formulation, employing equalorder bilinear ve...
متن کاملConditional simulations of flow in randomly heterogeneous porous media using a KL-based moment-equation approach
In this study, we extend the KLME approach, a moment-equation approach based on the Karhunen–Loève decomposition (KL), developed by Zhang and Lu [An efficient, higher-order perturbation approach for flow in randomly heterogeneous porous media via Karhunen–Loève decomposition. J Comput Phys 2004;194(2):773–94] to efficiently incorporate existing direct measurements of the log hydraulic conductiv...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2003